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Abstract 

In the presence of a non-vanishing chemical potential the eigenvalues of the Dirac operator be- 
come complex. We calculate spectral correlation functions of complex eigenvalues using a random 
matrix model approach. Our results apply to non-Hermitian Dirac operators in three-dimensional 
QCD with broken flavor symmetry and in four-dimensional QCD in the bulk of the spectrum. 
The derivation follows earlier results of Fyodorov, Khoruzhenko and Sommers for complex spectra 
exploiting the existence of orthogonal polynomials in the complex plane. Explicit analytic expres- 
sions are given for all microscopic fc-point correlation functions in the presence of an arbitrary even 
number of massive quarks, both in the limit of strong and weak non-Hermiticity. In the latter case 
the parameter governing the non-Hermiticity of the Dirac matrices is identified with the infiuence 
of the chemical potential. 



1 Introduction 



Random Matrix Theory (RMT) does not only provide a very useful tool for spectra of real variables 
such as energy levels. Also in the case where the eigenvalues of the underlying Hamiltonian become 
complex RMT describes universal properties in a variety of different physical models such as localiza- 
tion in superconductors [||], dissipation and scattering in Quantum Chaos ^, ^ |5| or chiral symmetry 
breaking in Quantum Chromodynamics (QCD) |^]. It is the latter subject which has motivated the 
results presented here. In the presence of a non-vanishing chemical potential for the quarks the eigen- 
values of the Dirac operator become complex. The regime with broken chiral (or flavor) symmetry 
which we want to describe is essentially non-perturbative in nature. One of the prominent techniques 
starting from first principles is to perform Monte-Carlo simulations on a lattice. However, in the 
presence of complex eigenvalues the simulations run into serious problems due to the presence of a 
complex determinant of the Dirac operator, which have not yet been overcome to large extent (for a 
review see and references therein). A better analytical understanding of the correlation functions 
of Dirac eigenvalues in the vicinity of the origin which are very sensitive to the symmetry breaking 
would therefore be very useful. The question of having better control over the lattice simulations in 
a certain regime is not of academic nature since the chiral phase transitions is now also investigated 
experimentally. 

In the past few years RMT has been developed as a powerful analytic tool to investigate this topic 
(for a recent review see (P). The reason for its applicability has by now been well understood on 
field theoretic grounds, analytic expressions have been provided for correlation functions at vanishing 
and non-vanishing temperature and even the qualitative picture of the QCD phase diagram has been 
understood in terms of RMT. However, analytic expressions for microscopic correlations in the presence 
of a chemical potential have not been derived so far and the present article is meant to partially close 
this gap. The only exception is the quenched limit with zero quark flavors. We will recall that the 
quenched correlations follow from ^ for weakly and from |10| for strongly non-Hermitian spectra. The 
nearest neighbor spacing distribution is also known in both cases, and respectively, but not in a 
closed form. It has been already compared to quenched QCD lattice data with chemical potential in 
|11] and a cross-over has been observed. Finally there are results for the spectral rigidity following 
semi-classical arguments. 

We will focus on the most simple model, the Hermitian matrix model or Unitary Ensemble, and 
discuss its generalization to incorporate complex eigenvalues in the presence of massive quarks. This 
model is relevant for Euclidean QCD in three dimensions (QCD3) jl^ and for the correlations in the 
bulk in four dimensions (QCD4) as observed in |0, 15|. Because of the absence of chiral symmetry in 
odd dimensions the corresponding global symmetry is flavor symmetry, which can be spontaneously 
broken from U{2Nf) U{Nf) x U{Nf) for an even number of quark flavors 2Nf. Recently evidence 
has been provided from lattice data |]l6| for the existence of a non- vanishing condensate in QCD3. In 
the following we will assume to be in the broken phase. In particular we will not address the question 
of different phases at high density from the RMT point of view []l7| . 

The corresponding correlation functions at vanishing chemical potential have been first derived 
from RMT for massless [13] and massive fiavors |18] and proven to be universal It has been 

observed that they can be related to finite- volume partition functions |2C]. Very recently all correlation 
functions have then been derived entirely from the underlying chiral Lagrangian using supersymmetry 
|2l|| and using the replica method |^], including the corresponding sum-rules |]22| , p3| . These results 
justify the RMT approach a posteriori and put it on firm field-theoretic grounds. 

Our approach will be more phenomenological here as it does not yet follow from a chiral Lagrangian. 
In complete analogy to QCD4 a chemical potential Hchem can be introduced on the level of the QCD3 
Lagrangian for the quarks, leading to a shift of the time derivative inside the Dirac operator do — > 
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do + fJ-chem- This Step clearly renders the Dirac eigenvalues to be complex. In the RMT formulation 
of Euclidean QCD4 at Hchem = the Dirac operator is replaced by an off-diagonal block matrix of 
complex matrices, offdiag(H^, W^). Its eigenvalues are real and occur in pairs of opposite sign, reflecting 
chiral symmetry. The insertion of a chemical potential is then mimicked by shifting both blocks by 
p^(t) _> w^^^ + ifJ-chem- After this step, the RMT eigenvalues are complex. So far no analytic results 
for the microscopic correlation functions have been obtained as a function of fJ-chem- 

In order to make some progress we address the RMT of QCD3 which is much simpler. Here, the 
Euclidean Dirac operator is replaced by a random matrix H which is Hermitian. The flavor symmetry 
breaking becomes visible after switching on quark masses (for a detailed discussion of symmetry 
breaking in QCD3 see e.g. In QCD3 these masses are purely imaginary and have to occur 

in pairs of opposite sign [p^], the quark mass matrix reading diag(zmi, . . . , irriMf, —imi, . . . , —irriN^). 

The full determinant of the Dirac operator is therefore given by det[-ff — inif] det[-fr + imj] in 

the presence of 2A^j quark flavors. As a first attempt we introduce the chemical potential ^chem i'^ 

the 

same way as it was proposed in the RMT for QCD4 by simply shifting H H + i^chem- However, this 
shift is trivial^ since the determinant can still be diagonalized by a unitary transformation. The Dirac 
operator eigenvalues simply obtained a constant imaginary contribution, xj — > Xj + ificherm where the 
Xj are the original real eigenvalues. This can also be seen from where a more general model with 
H + iV = J was studied for Nf = 0, with T being a diagonal matrix. The imaginary part of J was 
fixed with a delta- function constraint in the measure. Specifying F = ^chemX immediately leads to a 
constant shift of the real eigenvalues into the complex plane. The equivalence of the delta-measure 
and the canonical exponential measure was shown in |24| in the microscopic limit. 

When we take the large-A^ limit to determine the RMT correlations there are two possibilities. 
If we keep N Hchem fixed similar to the microscopic rescaling of eigenvalues and masses, Nxj and 
Nnif respectively, we reobtain the known massive correlations at shifted eigenvalues. The second 
possibility is to keep fJ-chem fixed at ^ 00 which immediately leads to a complete quenching and 
thus to the Nf = correlations [^] . Although in general a shift of -ff by a constant matrix known as 
random plus deterministic correlations may drastically change the spectrum (e.g. this is not the 
case here introducing Hchem iii the most naive way. 

In order to incorporate the influence of a chemical potential that leads to complex eigenvalues we 
will thus pursue another way. Namely we will replace the Hermitian matrix H hy a complex matrix 
J, without changing further the structure of the Dirac operator determinant. Both Hermitian and 
anti-Hermitian part of the matrix J will be averaged with a Gaussian measure. Given this general 
structure we are left with a few questions. Following there exist two different large-A limits, the 
limit of weak and strong non-Hermiticity. Both limits provide different correlation functions already 
in the flavorless case as discussed in In the following we will generalize these results in both limits 
to an arbitrary number of massive quark flavors in sections ^ and ^ It is clearly the first limit of weak 
non-Hermiticity which makes contact to the problem posed. In this limit there exists a parameter 
a which interpolates between non-Hermitian (a > 0) and Hermitian matrices (a = 0). The role of 
a as a function of ^chem will be discussed at least qualitatively. The microscopic large-A limit we 
will take slightly differs from Q in the weak non-Hermitian case. Because of the analogue []l3| of 
the Banks-Casher relation relating the condensate to the Dirac operator spectral density (per unit 
length) at the origin we are interested in the eigenvalues very close to zero. We will therefore define a 
microscopic origin scaling limit in which both the real and imaginary part of the complex eigenvalues 
Zj of J will be equally rescaled by the mean level spacing. In the limit a — > we will then recover the 
correlation functions of the Hermitian model |13, 18 1 with 2Aj massless or massive flavors. 

^Similar to shifting by a real constant c, H H + c, which leaves the bulk correlations (for Nf = 0) unchanged 
due to translational invariance together with the universality of the measure. 
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The article is organized as follows. In section |2| we present the general method how to calculate 
complex eigenvalue spectra with massive quarks. Here, we follow closely Q using orthogonal polyno- 
mials in the complex plane and show how to extend their results to Nf ^ 0. In the sequel we then 
discuss the two microscopic large- A'^ limits of weak and strong non-Hermiticity in sections ^ and ^, 
respectively. Section ^ contains a discussion of the universality of our results and a comparison to the 
Hermitian case. We conclude in section |5|. 

2 General result: massive correlation functions at finite- 



We start this section by stating the RMT partition function for QCD3 with a complex Dirac operator 
and 2Nf massive flavors 



■'QCD3 



{{ruf}) 



Nf 



dJdJ^ Yl \det[J -imf]f exp 
/=i 



N 



Tr ( JJ^ - T^eJ^ 



(2.1) 



Here, J is a complex matrix and we integrate over independent matrix elements. We have chosen to 
take the absolute value of the Dirac operator determinant in order to make the partition function real 
valued. It can be written in terms of complex eigenvalues 



4'^^)({m^}) = X{[d%X{\z,-imf\^ w\z,)\\^n{zi^...,zn)\\ 

j=i V /=i / 



(2.2) 



with the Vandermonde reading An{zi, . . . , zm) = I\k>i{zk — zi). The weight function is given by 

N 



w'{z) 



exp 



^{z' + z 



(2.3) 



The diagonalization of J has been described in detail in |25] for Nf = and we only repeat here 
the ingredients important for the Hermiticity properties. To justify the form of the average in the 
partition function eq. ( |2.1| ) we start from the following decomposition of an arbitrary complex matrix 



J 



'1 



1 + T 



-A 



(2.4) 



Here, H and A are both Hermitian matrices with Gaussian weight and equal variance (1 + t)/{2N). 
The parameter r G [0, 1] governs the degree of non-Hermiticity. While for r = we have maximal 
non-Hermiticity for r = 1 we obtain back a Hermitian matrix J = H. Furthermore, r controls the 
correlations between matrix elements of J as < JkiJik > = t/N and < JkiJki > = At large-A^ 

we will either keep 1 — r or N{1 — t) fixed, the strong or weak non-Hermitian limit. The product of the 
Gaussian measures 'P{H) ~ exp[— A^/(l -|- r)Tri?^] for the matrices H and A can be simply rewritten 
as a single measure for the complex matrix J, V{J) ~ exp[— A^/(l — r^)Tr(JJ^^ — rSReJ^)]. This leads 
to the appearance of particular the weight-function eq. ( |2.3|) for the complex eigenvalues Zj of J. 

The crucial observation is that the correlation functions for the partitions function eq. ( p.2| ) can 
be obtained using the powerful technique of orthogonal polynomials. Let us therefore define a set of 
polynomials p'^^^\z) orthonormal in the complex plane z = x -\- iyhy 

Nf 



j d\ n |z - im/|2 w\z) '^^)(z)pf ^^)(z* 
/=i 



^ki 



(2.5) 
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where cPz = dxdy and x, y € R. In the case Nf = these are given by appropriately rescaled standard 

Hermite polynomials, P^\z) ~ H}^{z^J N /t), an observation made in [^]. Furthermore we define the 
/c-point correlation function for finite- in the usual way: 

R^N^^\zi, ■■■,Zk) = J . . . (fzN VNiizj}) , 

1 T-r I T-r , -12 „.,2 




'PnHzj}) = [[{ [[\zi-imf\' w'{zi)]\AN{zi,...,ZNr . (2.6) 



Using the standard technique of orthogonal polynomials |27| the correlation functions can be expressed 
in terms of the kernel of the orthogonal polynomials 

Nf 7V-1 
i^^'^^)(zi,4) = n \zi-imf\\z*,+imf\ w{z,)w{zl) ^ p/'^^)(zi)P/'^^^(4) (2-7) 



as 



/=i 1=0 



R^P\zu...,Zk) = det \k'^'''\z„z*)] . (2.8) 

i,j=l,...,k L ■'J 

In order to calculate the correlation functions for an arbitrary number of massive flavors 2Nf we would 
have to determine the orthogonal polynomials in eq. ( |2.5D and then evaluate the kernel. This could 
in principle be done in analogy to [l^ starting from Hermite polynomials for Nf = and iteratively 
adding massive flavors (see also Theorem 2.5 in |^^). We will proceed in another way by directly 



relating the massive correlation functions R^^^^^ to the flavorless ones, /?yy^ , SiS it was done in [g9[ for 
the chiral ensembles of real eigenvalues. The advantage of this derivation is twofold. First of all we will 
immediately obtain a closed and simple form for all Nf, without calculating intermediate polynomials. 
Second, the RMT universality of our results will be directly inherited given the universality of the 
flavorless case can be shown. 

We now describe how to obtain the massive correlation functions. The key observation borrowed 
from 1 29] is to incorporate the Dirac operator determinant with 2Nf flavors into a larger Vandermonde 

[[[[\zj-imf\ \An{zi,...,zn)\ = ' -—^2 —■ (2-9) 

j^if^i \^Nf\imi, . . . jiiriNf)] 

Note that due to the pairing of imaginary masses in QCD3 we do not need to impose an extra degen- 
eracy to the power of the Dyson index /3 as it was done in |29] for the chiral ensembles. Furthermore, 
the relation ( |2.9D is the technical reason why we have to take the absolute value of the complex Dirac 
operator determinant in eq. (2.1) to begin with. Otherwise we would not be able to relate the flavorless 
correlations to the massive ones. Inserting the identity (|2.9| ) into the definition (^|^) we obtain^ 

-r(0) Nf 0(0) / • • \ 

„i2Nf), _ N\ ^N+Nf U 1 R)^'+N^{zi, . . . , Zk,imi, . . . ,imNf) 

" - (iV + iV/)!42^/)({^^})jL\u;2(,^^) |A^,(zmi,...,im^,)P 

(2.10) 

In order to evaluate the mass-dependent normalization factor ^j^^^^({m^}) we also insert eq. ( |2.9D 
into the definition (p. 21) to obtain 



^ ~ {N + Nf)\f}^w^{imf)^^+^f\AMf{imi,...,imNf)\''- ^ ' 



Strictly speaking eqs. (|2.10[)-(|2.12|) only hold for finite-A'' in the normalization where the weight function eq. 



(2.3) is A''- independent (see ESfl). In the large- A'' limit to be taken later this difference becomes immaterial. 
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With these two results we obtain the following master formula for the massive fc-point correlation 
function of complex eigenvalues in terms of the known flavorless correlation functions of [^] : 

(2^) R'-°l^{zi,...,Zk,imu...,imNf) 

Rn {zi,...,Zk) = ^^oy- ^ . (2.12) 

Ry^,{imi, . . . ,imNf) 

The right hand side is entirely given in terms of the zero-flavor kernel eq. (^]^) from which we 
explicitly display here in terms of the Hermite polynomials, 

kP{zi,z*,) = w{z^)w{zl)—=== V —Hi I \r-z^ I Hi I I • (2.13) 




Since all arguments on the right hand side of eq. ( p. 12 ) are complex the insertion of the imaginary 



quark masses imf does not constitute an analytical continuation, in contrast to the case of real 
eigenvalues considered in |^^. We thus do not have to face the subtleties that occured there when 
taking the absolute value in the large- limit. 

Let us add that we can use exactly the same strategy when restricting the partition function eq. 
( p.2| ) to real eigenvalues. We thus recover without any effort from eq. ( p. 12 ) the results of [18| for 



massive microscopic correlation functions out of the flavorless ones. In fact we will reobtain their 



results from our correlation functions of complex eigenvalues in the Hermitian limit in subsection 3.1 
While eq. ( ^.12| ) establishes an exact result for finiteS-A'^ we will be primarily interested in taking 
the microscopic large- limit. As it has been shown Q there exist two distinct microscopic large- A^ 



limits: the limit of weak and strong non-Hermiticity. The explicit form of the correlations eq. ( 2.12 ) 
in these two limits will be the subject of the next two sections. 

3 Microscopic limit I: weak non-Hermiticity 

Let us start by by defining the limit of weak non-Hermiticity recalling the results of |9|, We will 



take the limit r ^ 1 in a controlled way in which the matrix J in eq. (2^) becomes almost Hermitian. 
Namely we will keep 

lim 2A^(1 - r) = (3.1) 

fixed in the limit where A^ — > oo. The Hermitian limit can then be recovered by taking q ^ 0. From 
the requirement of keeping the weight in eq. ( p. 3D finite we immediately see that we have to rescale 

the arguments zi_2 of the kernel K^^^^\zi, Z2) in eq. ( |2.7D appropriately. The authors of find 
non-trivial correlations when 

^e{zi - Z2) = 0{l/N) and 9m zi,2 = 0{1/N) , (3.2) 

with 1/A^ being proportional to the mean level spacing while ?R.e zi^2 = C(l) is kept finite. In other 
words they require SRe z ^ 9m z in addition to the limit eq. (p.l|). Let us furthermore define the 
mean spectral density of the real part of the eigenvalues only: 

u{x) = ^ < p{x = ^e z)> . (3.3) 

Here, the large-A'^ average is taken without rescaling the eigenvalues or unfolding, leading to the 
macroscopic density. For a Gaussian weight this is just the semi-circle z^(x) = ^V'i^ — ■, since the 



average over the matrix A drops out. In |25| the final result for the asymptotic kernel only depends 
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on the average real part of the eigenvalues, X 
apart from a pre- factor. It reads 



3f?e(zi + Z2)/2^ through the mean eigenvalue density 



K 



AT- 



N 



vra 



exp 

■KV{X) 
-■KU{X) 



-XiV9m(zi - Z2] 



exp 



— - y-^jia z\ + xsm. Z2 ^ 



du 



exp 



+ iN{zi - zl)u 



(3.4) 



After unfolding the kernel, the correlation functions and the eigenvalues according to eq. ( |3.2| ) the 
zero- and the A^y-flavor correlators can be immediately read off from eqs. (2.8) and ( p. 12 ), respectively. 

Because of the Banks-Casher relation between the condensate and the mean spectral density at 
the origin from QCD3 [13| we are interested in the microscopic origin scaling limit. We can therefore 
relax the condition 5Re z ^ Q'm z since we rescale both real and imaginary part in the same way. The 
microscopic origin scaling limit is defined by introducing variables 



N{^(i z + i9m z) 



Nz =^ 



Nmt 



(3.5) 



which is kept fixed at ^ 00. In addition we have rescaled the masses accordingly since they appear 
on the same footing in the determinant of the partition function eq. (^]^) . In consequence the kernel 
eq. ( |3.4D simplifies furthermore as now the average X is also quantity of order 0{1/N) vanishing at 
large- A^. We thus obtain for the zero-flavor kernel of weakly non-Hermitian eigenvalues 

1 



4°^ (6,^2) 



7ra 



■ exp 



where0 



7rj/(0) 



du 



exp 



(3.6) 



(3.7) 



Note that ga{C) is an even function in ^ due to the range of integration. Although the result eq. (|3.4D 
was obtained with a Gaussian potential only, with 7rz^(0) = 1, we keep the dependence on J^(0). It 
directly connects to the condensate at zero chemical potential and in the Hermitian limit it plays the 
role of a universal parameter. The issue of universality will be discussed separately below. 

The result for the zero-flavor A:-point correlation functions simplifies as the exponential pre-factors 
of the kernel can be taken out of the determinant. Defining the rescaled or unfolded correlation 
functions, 

pf{ii,...,ik) ^ lim N-^R%\N-^i^,...,N-\u) (3.8) 



we obtain from inserting eq. (|3.6D into eq. d2 



1 

vra 



k k 
]Jexp 

i=l 



2 

a 



l<],l<k 



Ci)] 



(3.9) 



where we have used that Q'm^^ = 9m^^*. This also gives us the building blocks for the denominator 
in eq. ( p.l2| ) with purely imaginary arguments. We can now read off the massive correlation function 
using the zero-flavor kernel eq. ( ^.6| ) in eq. ( |2.12| ) 



Prk6,...,e.) 



det 






-iPh) 


l<j,l<k 








l<f,h<Nf 






-iph) . 



deti 



i<f,h<Nj[Ks\ilJ'f, -iph)] 



(3.10) 



Note the difTerence by a factor of i in the definition compared to 
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Figure 1: The flavorless density in the complex plane pi {^) at = 0.2 (left) and the microscopic 
density of real eigenvalues for 2Nf = 2 massless flavors (right). 



The first main result of this article, the massive A'^y-flavor fc-point correlation function is thus reading 



Ps (6,...,Cfc) 



vra 



k k 



^ n 



exp 



1=1 



-^9m 



det 


9a{ij - 


CD 


gaiCj +iph) 


l<i,l<k 








i< .f,h<Nf 


_ ga{ipf 


-er) 


ga{ipf + iph) _ 



deti<f^h<Nf[gaiifJ'f + iph)] 



(3.11) 

We mention as an aside that from the connected part of the two-point function we can determine the 
absolute value squared of the corresponding massive kernel, |-f^5^^''^('?i, Cl)!^- This follows from the 



definition {2.1 



In the following we will give a few explicit examples and discuss the limit of massless flavors. For 
Nf = 0, 1 the microscopic density reads 



1 



ira 



exp 



p's'iO = — exp 
vra 



o 2/- 
— ^^m t, 



5„(2z9=m i) , 



ga{2iQm ^) 



ga{C + iP')ga{ip- C) 
ga{2ip) 



(3.12) 
(3.13) 



While the flavorless result eq. (3.12) follows from the origin scaling limit of |3^ the Nj = 1 flavor case 
is new. The massless limit is easily obtained by setting /x = in eq. (|3l^ ). We have plotted in Figs. 
I and I both cases Nf = and 1, respectively, for a comparison to the known Hermitian limit. 

With the general structure of the massive correlators being clear from eq. ( ^.11 ) we also give 
explicitly the massless two-flavor case: 



gQ,(2i9m ^) 



ff«(g)ga(r) 

5a (0) 



+ 



(yie-'/2sin^-e5a(0) (yfe-"V2sine-r5a(r) 



(3.14) 



8 



In contrast to the correlation functions of real eigenvalues |13] the massless correlation functions do not 



simplify here. The procedure of sending successively all masses to zero in eq. ( p. 11 ) leading to higher 



order derivatives of Qa will not simplify the determinant structure. The reason for this will become 
clear when comparing to the Hermitian limit below. Taking for example the microscopic density with 
Nf massless flavors we will obtain a determinant of size Nf + 1 even in the limit a ^ 0. It is only 
through non-trivial identities that this expression can be rewritten as a sum of bilinear combinations 
of half-integer Bessel-functions for an arbitrary Number 2Nf of massless flavors as first derived in |jl^. 

3.1 The Hermitian limit 

We will now discuss the Hermitian limit when a ^ 0. In the simplest case of zero flavors the 
microscopic density of the Hermitian matrix model is simply constant. For the complex density eq. 
( ^.12 ) this feature is maintained on the real axis while the density smoothly spreads into the complex 



plane as depicted in Fig. || (left) (see also Taking the limit a ^ of eq. ( p.l2| ) we obtain a 

delta- function times a constant, the density of the real eigenvalues: 

limpl°^(0 = ^5 fV29m ^) lim c/c.(2i9m = z^(0)(5 (9m ^ , (3.15) 

a— >0 a/vT V / a^O 



where in the last step we have used lima^o 9a (0) = v27ri^(0). If we unfold the microscopic density as 
ps{C) = liiHAT^oo -Di?7v(-C'^) with respect to the full mean level spacing, D = l/(7ri/(0)A^), we obtain 
the parameter free result Pg\^) = 1/vr [13]. Next, we turn to 2Nf = 2 where we take the massless 



limit for simplicity. In the Hermitian model the density is given by [13| 



pL''(f) = . i^l^, (3.16) 




(2) 

Figure 2: The microscopic density in the complex plane ps (0 of 2Nf = 2 massless flavors for 
different values of o? = 0.5 (left) and = 0.1 (right). 



which is plotted in Fig. JT| (right). The corresponding density in the complex plane, eq. ( 3.13| ) with 



= 0, is plotted in Fig. |^ for two different values of a. It very nicely combines both properties of Fig. 
|l|, the oscillations of eq. ( |3.16| ) on the real axis (right) and the same spreading into the complex plane 
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(left) of eq. ( 3.12 ) for Nf = 0. For decreasing a the distribution gets narrower, slowly approaching 
the delta-function. This is also the reason for the different normalization in Figs. |^ and ^ for different 
values of a. More mathematically we obtain from eq. ( |3.13| ) 



z^(0)5(9m 1 



sin2(7ri/(0)O\ 
(MOW ) 



(3.17) 



When unfolding with respect to the mean level spacing D we thus recover eq. ( 3.16 ). 

Let us now turn to the Hermitian limit in the general case. For that purpose we first consider the 
Hermitian limit of the zero-flavor Kernel eq. (|3.6| ) as a building block. We obtain |^ 



hmXf (^1,^2*) 



iO)V2 6 



Qm^^i + 9m2^2 



sin(^z.(0)(6-C2*)) 

Mo)(ei-e2*) 



(3.18) 



After unfolding with respect to D instead of 1 /N this becomes precisely the parameter free universal 
sine-kernel of ||3^ on the real axis. It therefore does not come as a surprise that we also recover the 
known result for the microscopic massive fc-point correlations 



limpr'^(6,...,e.) 



(3.19) 



1=1 



det 

l<j,l<k 
l<f,h<Nf 


" sin(7ri/(0)(^j - ^0) sm{7r ly (0) {^j + i^h)) 
sin(7ri^(0)(i^/ - Ci)) sinh(7ri/(0)(^/ + Hh)) 




deti<f^h<Nf 


sinh(7rzy(0)(/i/ + l^h)) 





Taken in units of D it exactly coincides with the corresponding result for the Hermitian model as 



derived in |21] from the underlying field theory. The equivalence with the original form |18] can be 
shown using consistency conditions on QCD3 finite volume partition functions 

Let us come back to the remark at the end of section |2|. Had we taken the partition function eq. 
( p.2| ) for real eigenvalues in the beginning, together with a general polynomial in the exponential of the 
weight function eq. ( |2.3[ ), we would have immediately arrived at eq. ( 3.1S| ) after employing eq. ( 2.12 ) 
and the zero-flavor kernel eq. (3.18). We therefore rederive the results of |l^ without any effort. The 



universality of our results then follows from that of the sine- kernel |32, [19[ . 
3.2 Universality 

Up to now we have only discussed a large-A'^ limit in which the eigenvalues are rescaled with the 
mean level spacing, the microscopic limit. When we want to discuss the issue of universality we have 
to distinguish another type of large-A^ limit, the macroscopic limit. In this limit the short range 
fluctuations typically of the order 1/A^ are smoothed and the macroscopic (or smoothed or wide) 
correlation functions are obtained. Technically speaking the large- limit is taken without a rescaling 
(or unfolding) of the eigenvalues. This limit is also meaningful since here universal correlation functions 



can be obtained |33] as well. An example for such a macroscopic correlation function is the density in 
eq. (|3.3|), the smoothed eigenvalue density of the real part of the eigenvalues. For all three Gaussian 
ensembles of real eigenvalues it is given by the Wigner semi-circle. However, for a measure including 
terms of higher order than Gaussian the macroscopic spectral density z^(x) itself is non-universal as it 
explicitly includes all parameters of the higher order terms. 
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Going back to complex correlation functions the authors of |25] have shown a certain degree of 
universality for the macroscopic spectral density p{x,y), with z = x + iy, in the limit of weak non- 
Hermiticity. Using supersymmetry they have shown that it only depends on the real part x trough the 
combination z^(x), eq. (|3.3|). Hence the macroscopic density as a function of y only, denoted by Px{y) 



in [25 1, is universal in the sense that when the scale is set by fixing the real part of the eigenvalues as 
X = X in [25 1 or x = as in our case of origin scaling, plays the role of a universal parameter 
in the case of a non-Gaussian weight function. Similar results have been obtained in [Q| for weakly 
non-symmetric real matrices. 

A more interesting point is the universality of the fluctuations in the microscopic correlations we 
have calculated in eq. ( ^.11 ). We have shown that given the flavorless correlations determined in 
B are universal their universality would automatically carry over to our massive results trough eq. 



( ^.12 ) in the microscopic limit. From the derivation [^5| of the kernel eq. (|3.4| ) it is obvious, that the 
occurrence of the mean spectral density of the real eigenvalues inside the integral eq. ( |3.4| ) is 

strictly related to the form of the measure eq. . We can therefore only conjecture that the kernels 
eq. ( |3.4| ) and (^]^) and thus the corresponding correlations are universal also for a more general weight 
than eq. (^^). However, we can offer a non-trivial check in the Hermitian limit a ^ 0. As we have 
seen in this limit the flavorless kernel eq. ( |3.4[ ) maps to the sine-kernel eq. ( |3.18| ) which is guaranteed 
to be universal (3^, In the microscopic origin scaling limit the universal parameter z^(0) occurs 
in the correct place after taking the limit a ^ of eq. (^|^). The same argument thus holds for the 
Hermitian limit eq. ( 3.19| ) of eq. ( p. 11 ). It would be highly desirable to repeat the universality proof 
of [^] for orthogonal polynomials in the complex plane in order to obtain eq. ( |3.(]| ) and consequently 
eq. ( |3.11| ) for an for an arbitrary weight function. 

For a very recent discussion of universality of complex spectra from the point of view of the 



Fokker-Planck equation we refer to [34 



3.3 Discussion of a as a function of fichem 

As it has been discussed already in the introduction the naive introduction of a chemical potential 
by writing the Dirac operator as -ff -|- ipchem does not lead to nontrivial correlations of complex 
eigenvalues. It either lead to a complete quenching or, after rescaling with keeping Nfichem fixed, to a 
constant shift in the eigenvalues. We therefore have introduced a truly complex Dirac operator J and 
we would like to understand the role of the parameter a governing the non-Hermiticity as a function 
of chemical potential Pchem- It would be therefore very instructive to compare to the RMT model in 
four dimensions [|6[ and identify the parameters. However, in QCD4 much less is known. Only the 
macroscopic density and its support are known [^ as a functions of Pchem in the RMT model while 
the microscopic correlations to be compared with lattice simulations are still lacking. 

What we can do is to calculate the variation of the imaginary part of the eigenvalues, 9m ^, in the 
limit a ^ and comparG it to ci vanishing potcnticil j^chem 

— > in [^. For simplicity let us consider 
the zero-flavor case with microscopic density eq. (|3.12| ). We can calculate the second moment along 
the imaginary axis < 9m^^ >= J^^dQm^ 9m^^ P's^ {0 given by 

< We > = {l + ^a^^ . (3.20) 

On the QCD4 side we can ask with which power of Hchem the unsealed imaginary partQ of the eigenvalues 
vanishes on the boundary of the support. It follows from that 

5m z = M ,\,^\/4-3f?e^z + . (3.21) 

*Note that in 0] for vanishing /ichem the Dirac eigenvalues are chosen to be purely imaginary. 



11 



We notice that the mean spectral density naturally occurs. It is therefore tempting to identify 
at least to lowest order 

aifJ^chem) ~ fJ'lhem ■ (3.22) 

Of course we do not have to make this identification. We could also simply take a as a fit parameter 
to compare with numerical data for the QCD3 Dirac operator with chemical potential once they are 
available. The analytical form of the microscopic correlations would still be fixed by eq. ( p. 11 ). 

Let us mention a direct application for QCD4 with chemical potential. It has been found empirically 
that also away from the origin the bulk correlations of Dirac operator eigenvalues can be very well 
described by RMT |l|] (for the most precise tests we refer to [|l5l). In the bulk of the spectrum the 
Dirac determinant and thus the chiral properties are no longer seen which makes our non-chiral RMT 
eq. (^]^) applicable. Thus we can take the zero-flavor microscopic correlation function of complex 
eigenvalues eq. (|3.9| ) to test QCD4 correlation functions with chemical potential in the bulk of the 
spectrum. As it was shown in the macroscopic spectral density in the complex plane is independent 
of the real part Ke z of the eigenvalues. This is consistent with the quenched microscopic density eq. 
( p. 121 ). One could thus check if the complex eigenvalues decay microscopically as depicted in Fig. Q 
(left) along a fixed line of of constant ?R.e z parallel to the real axis. Comparing for different values 
of Hchem would then fix the constant a as a function of fichem- We stress that also the higher order 
correlation functions are available analytically in eq. (3^). The corresponding real correlations eq. 
(f3.19) have been already tested in |15] for vanishing chemical potential. 



4 Microscopic limit II: strong non-Hermiticity 



In this section we treat the large- limit in which the parameter r, or more specifically the combination 



^(1 — t)/(1 -|- t) in the definition ( |2.4| ) is kept fixed. Here, there will be no limit possible to recover the 
correlations of the original Hermitian model and thus QCD3 with real eigenvalues, such as discussed 



in subsection 3.1. However, we will find nontrivial generalizations of the original work of Ginibre jlC 
It turns out, that when introducing quark fiavors the microscopic spectral density starts to decay 
compared to the zero-fiavor case, where it is entirely fiat on the support [|l^ (see eq. ([4.3| )). 

The major difference to the weak non-Hermitian limit comes from the fact that the mean level 
spacing is changed to be -D ~ 1/ \/]V. When looking at the weight function eq. ( |2.3| ) with 1 — r fixed 
we can define the microscopic origin scaling limit for strong non-Hermiticity as 



N(JR.e z + i9m z) 



(4.1) 



The microscopic kernel which is now defined as ^^5(^1,^2) = Kn{S^i/ ^/ N ,Q/ \/ N) / \/ N can be directly 
obtained for Nf = by taking the limit (4.1) of the result for finite- A'', eq. ( |2.13| ) p5[ |: 



1 



7r(l-T2) 



exp 



1 



2(1 - r2) V^^'^'* " " + 2 " + ^2 - ^2 



At equal arguments we thus obtain a constant microscopic density for zero flavors [10, 25] 



7r(l-T2 



■ (4.2) 



(4.3) 



It coincides with the smoothed or macroscopic spectral density of the unsealed eigenvalues |1C] which 
has been reobtained in using the supersymmetric method 



lim -^RMiz) 



1 



o 2 
xsva. z 



7r(l-r2) ' 
, 



^(?z 

' (l + r)2 + (1-t)2 
otherwise 



< 1 



(4.4) 
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The latter is independent of the number of quarks present as it is obtained from saddle point analysis 
where the quark determinants are subleading. We recall that in the RMT of QCD4 the support of 
the macroscopic density is as well given by an algebraic curve and, although not being flat, the density 
is independent of the real part of the Dirac operator eigenvalues. 

Next we turn to the microscopic correlation functions in the presence of 2Nf massive quarks. 
Inserting the zero-flavor kernel eq. ( |4.2| ) into eq. ( |3.1[1| ) now in the strong non-Hermitian limit eq. 
( |4.l| ) we obtain for the general fe-point function 



Ps 



7r(l 



k k 

J|exp 

1=1 



det 

l<j,l<k 
l<f,h<Nf 



exp 




exp 





exp 



exp 



det 



l<f,h<Nf 



exp 



1 



l-T- 
/^/ Ph 



1 



(4.5) 

Here, we have taken those parts in eq. ( |4.2| ) which factorize out of the determinant, leading to a partial 
cancelation. Let us explicitly display the simplest examples, the density with one massive flavor 



pS?'({) 



1 



7r(l 



1 



exp 



and the density with two massless flavors 



7r(l 



1 + 



exp 



(4.6) 



(4.7) 




Figure 3: The massless microscopic density in the complex plane ps ^ (^) for 2Nf = 2 (left) and 
27V/ = 4 (right), both at r = 1/2. 

In both cases it is no longer constant. For vanishing masses both vanish at the origin ^ = due to 
the level repulsion of the additional quark flavors. For large arguments they both approach the value 
of the macroscopic density eq. ( [4.41 ). As we can see in Fig. ^ for two flavors (right) the level repulsion 
is stronger. 
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We finally note that the microscopic correlations in the limit of strong non-Hermiticity can be 
obtained from those in the weak limit by taking a — > oo, as being mentioned already in |25]. In this 



limit eq. ([4.5| ) follows from eq. ( ^.11 ) when identifying = (1 — r"^). In the simplest example eq. 
( p. 12 ) leads to the constant density eq. (4.3). We thus find a crossover from the regime of weak to 
strong non-Hermiticity with growing a. A similar crossover has been seen for the nearest neighbor 
distribution in quenched QCD4 lattice data [O] as a function of growing fj^chem- 



5 Conclusions 

We have derived analytic expressions for all correlation functions of complex Dirac operator eigenvalues 
in the presence of an arbitrary even number of massive quark flavors. Here, we have replaced the 
Dirac operator for QCD3 by a general complex matrix as we have argued that a shift of an originally 
Hermitian operator by a constant chemical potential, ifJ-chem, would not lead to non-trivial correlations. 
We have solved the problem in the limit of weak and strong non-Hermiticity and have identified the 
parameter for non-Hermiticity a as modeling the influence of a chemical potential. Our results may 
also be useful for QCD lattice data with chemical potential in four dimensions in the bulk of the 
spectrum. In the limit of a Hermitian Dirac operator we could nicely reproduce the known results 
for real eigenvalues. Moreover, we could offer a very elegant alternative of calculating massive real 
eigenvalue correlations and proving their RMT universality. 

There are several open questions we have not addressed so far. The main task is a generalization 
of our results to chiral RMT as an effective model for QCD4 spectra close to the origin. The crucial 
step would be to find the corresponding orthogonal polynomials in the complex plane to make our 
techniques applicable. An open problem within our calculations is the question of universality of the 
correlation functions which we could only conjecture. The problem is to show that for an arbitrary 
weight function we will obtain orthogonal polynomials with the same asymptotic as the Hermite 
polynomials. Given the success for real eigenvalues in this respect progress should be possible. Another 
question is to find the finite-volume partition function corresponding to the partition function of 
complex eigenvalues. In other words one would have to find an appropriate unitary group integral 
given by the determinant of the massive kernel of complex eigenvalues. 

Another open problem would be the generalization to the other invariant random matrix ensem- 
bles, the orthogonal and symplectic ensemble, to incorporate complex eigenvalues. First results have 
been obtained only for the macroscopic spectral density of almost symmetric matrices using the su- 
persymmetric method. 
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